Random matrix model at nonzero chemical potentials with anomaly effects 
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Phase diagram of the chiral random matrix model with U(1)a breaking term is studied with 
the quark chemical potentials varied independently at zero temperature, by taking the chiral and 
meson condensates as the order parameters. Although, without the U(1)a breaking term, chiral 
transition of each flavor can happen separately responding to its chemical potential, the U(1)a 
breaking terms mix the chiral condensates and correlate the phase transitions. In the three flavor 
case, we find that there are mixings between the meson and chiral condensates due to the U(1)a 
anomaly, which makes the meson condensed phase more stable. Increasing the hypercharge chemical 
potential (^y) with the isospin and quark chemical potentials (/i/, ^q) kept small, we observe that 
the kaon condensed phase becomes the ground state and at the larger /iy the pion condensed phase 
appears unexpectedly, which is caused by the competition between the chiral restoration and the 
meson condensation. The similar happens when and are exchanged, and the kaon condensed 
phase becomes the ground state at larger below the full chiral restoration. 



CD 
(N 

> 

a\ 
o 
o 



X 



I. INTRODUCTION 

The phase diagram of QCD at finite density has long 
been challenged from the perspectives on the fundamen- 
tal properties of dense matter realized in the core of a 
compact star and relativistic heavy ion collisions. The 
chiral symmetry breaking in light quark sector is one of 
the crucial features in low-energy QCD, which should be 
reflected in any effective model. In the massless quark 
limit with flavor number Nf = 2 or 3, chiral symmetry 
breaking gives rise to continuously degenerate vacua and 
to emergence of the massless Nambu-Goldstone bosons 
as the most important degrees at low energies. In reality, 
nonzero u- and d-quark masses and heavier s-quark mass 
lift this degeneracy and explain the observed pseudo- 
scalar meson spectrum. 

A system at finite quark number densities is character- 
ized with the quark number chemical potentials /z/ (/ = 
u, d, s) in the grand canonical description. The simplest 
situation at finite density may be the case with equal 
chemical potentials ^ = /in = j-id ~ Ms for all flavors, 
where the s-quark density may be suppressed compared 
with the u- and d-quarks due to the mass difference. 
However, the quark chemical potentials yU/ are generally 
unequal in various physical situations. There are more 
protons than neutrons in nuclei and also in neutron stars 
because of the electric charge of protons. In the core of 
compact stars in beta equilibrium, the heavy mass of the 
s-quark may be overcome by the large electric chemical 
potential, which may result in the appearance of s-quark 
degrees of freedom. In ultra-relativistic heavy-ion colli- 
sions the total strangeness number is constrained to zero. 

The phase structure at finite chemical potentials has 
very rich theoretical possibilities Historically, in 

dense nuclear matter, nucleon superfluidity, pion and 
kaon condensations, and hyperon mixture were studied 



based on empirical nucleon interactions and chiral per- 
turbation theory Q. More recently, color superconductiv- 
ity was extensively investigated and the kaon condensed 
phase acquired renewed interests there Q. Possibility of 
inhomogeneous phases is also discussed recentljQ 01 • 

Another motivation to consider the phase diagram 
with unequal chemical potentials is that one can directly 
examine in lattice QCD simulations the Nf = 2 case 
with finite isospin and zero quark chemical potentials, 
/i/ = {jjLu - Md)/2 ^ and /x, = (m„ -I- /Xd)/2 = 0, where 
the measure for the importance sampling is real. We note 
that it may be possible to simulate the Nf = 3 case with 
/i/ 7^ and = = 0, at least in principle. 

The ground state at finite isospin chemical potential 
but at zero quark number density was studied with the 
chiral Lagrangian, which revealed that the pion conden- 
sate appears once /x/ exceeds half the pion mass m,r/2 
Q. Indeed, in the chiral limit, the degenerate ground 
state is completely chiral-rotated to the pion condensed 
state at infinitesimal external field /z/ 7^ 0. Adding a 
finite quark mass 7^ 0, which makes the pion mas- 
sive, we have a competition between the two alignments; 
(qq) and (U750?) directions for and fij, respectively. In 
the three-fiavor case, kaon condensation appears once the 
hypercharge chemical potential fiy exceeds the threshold 
given by rriK with /i/ =00- 

The chiral random matrix (ChRM) model is one of the 
models which share the chiral symmetry with QCD. Thus 
it will serve as a useful model providing qualitative fea- 
tures of the QCD phase diagram. The aim of this paper 
is to explore the phase structure of the (ChRM) model, 
including the possibilities of the meson condensations, as 



^ The p-wave pion condensed phase too accompanies 
inhomogeneity 0] . 
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a function of the quark chemical potentials fif (/=u, d, 
s) at zero temperature. 

In the previous work of the ChRM model with two 
flavors done in Ref. Q it was shown that the pion con- 
densed phase appears at finite fii above a critical value 
file as suggested by [1], whereas it eventually disappears 
as fii is further increased, which is understood as the chi- 
ral restoration at high density. It was also indicated that 
with small nonzero fij the chiral restoration occurs in two 
steps: as fiq is increased with (e.g.) negative < fixed, 
the (Jd) condensate discontinuously decreases first, and 
then the (uu) at slightly larger fiq. The two- (or three-) 
step transition along the fiq axis with finite /i/ and/or 
IJ,Y is also observed in the three flavor case of the ChRM 
model i, [3 

and in the two and three flavor cases of the 
Nambu-Jona-Lasinio (NJL) model [TT| - [l^ . 

This two-step transition happens evidently because 
these models do not include flavor mixing interactions 
and then each quark sector independently responds to 
each quark chemical potential fif. The U(1)a anomaly 
in QCD is known to induce the flavor mixing and tends to 
unify the two chiral transitions into one. This is demon- 
strated in the NJL model with U(1)a breaking effective 
interactions, which results in the usual phase diagram 
with a single transition[ll|. 

The U(1)a anomaly effect is indispensable to build 
a low-energy effective model of hadrons in medium as 
well as in vacuum (l6l - [l8j . Although in a conventional 
ChRM model the U(1)a anomaly effect is not treated 
appropriately so as to affect the nature of the chiral 
transition, recently we have succeeded in constructing 
the model in which the U(1)a breaking term generates 
the flavor- number dependence of the transition order [l9j. 
This model in the massless quark limit shows the second- 
order phase transition for Nf = 2, whereas the flrst-order 
one for Nf = 3. We explored the phase diagram of the 
model in the T-fiq-mq space with {qfQf) as the order 
parameter [20|. 

It seems now quite intriguing to investigate the phase 
diagram allowing the pion and kaon condensates in the 
ChRM model with Nf =2 and 3 as a function of the 
chemical potentials fif, with emphasis on the role of 
the U(1)a anomaly. Similar studies were previously per- 
formed employing the NJL model with U(1)a anomaly 
effects [11, [H, [lil. In Ref. Q the phase structure 
for Nf = 2 is studied only in small fij region below 
the threshold for the pion condensation, and the study 
is extended by including the pion condensed phase in 
Ref. [l^ . Ref. [2l| deals with the Nf = 3 case restricted 
at fij ^ and fiy = 0, where the possibility of kaon 
condensation is ignored. 

This paper is organized as follows. After a review of 
the QCD phase structure with finite chemical potentials 
in Sec. 2, we introduce the ChRM model with anomaly 
effects in Sec. 3. The model phase structure including 
the possibility of the mesonic condensation is investigated 
analytically and numerically in Sec. 4 and 5 for Nf = 2 
and 3, respectively. Sec. 6 is devoted to a summary. 



II. QCD WITH QUARK CHEMICAL 
POTENTIALS 

In this section we briefly review the symmetry break- 
ing pattern of QCD with nonzero quark chemical poten- 
tials. The QCD partition function with Nf quark flavors 
at finite temperature and chemical potential is regarded 
as the average of the Dirac operators with non-Abelian 
gauge field action 5ym: 

Zqcd= VAY[detiD{fif)+mf)e-^^^, (1) 
/ 

where D{fif) = "^^{d^ — igA^ + B^fif) is the (Euclidean) 
Dirac operator for a fiavor / and is the mass. The 
quark number chemical potential /i/ is introduced inde- 
pendently for each quark flavor / with a constant Lorentz 
vector = (0,0,0,1). 

In the massless limit ruf = with zero chemical po- 
tential fif = 0, the classical QCD Lagrangian possesses 
U(A^/)lxU(A^/)r chiral symmetry, which is broken down 
to U(l)BxSU(A^/)LxSU(iV/)R at the quantum level due 
to the axial anomaly. The U(1)b invariance results in the 
quark number conservation. Non-perturbativc QCD dy- 
namics is believed to break this symmetry spontaneously 
to U(1)bxSU(A^/)v, generating a large part of the nu- 
cleon mass and the light pious as the Nambu-Goldstone 
(NG) boson. In reality, nonzero quark mass term acts 
as an external field to select the ordinary vacuum state 
with nonzero (qq) condensate out of the nearly degener- 
ate vacua. 

The chiral symmetry of QCD is intact at finite quark 
number chemical potential as far as it is flavor-singict, 
which is usually assumed in the context of flnitc density 
QCD. Once the fiavor-nonsinglct potential is set nonzero, 
it gives a stress on the fiavor symmetric ground state in 
addition to the quark mass term. Symmetry breaking 
pattern in this situation has been studied most conve- 
niently with the chiral Lagrangianj^ Q. 

Let us consider the Nf = 2 case with the small u- 
and d-quark masses iriu = md- At zero chemical po- 
tential there appear nearly degenerate vacua connected 
with each other by the SU(2)xSU(2)~0(4) transforma- 
tion, and the true ground state aligns to the direction of 
the quark mass term with SU(2)v invariance. When the 
isospin chemical potential fii ~ (/i„ ^ fLd)/1 is applied, 
this invariance explicitly breaks down to U(l) £ SU(2)v 
around the 3rd isospin axis. The mass m and the poten- 
tial fii compete to fix the ground state. For fij > TO,r/2, 
a new class of the degenerate vacua appears with nonzero 
pion field and the ground state spontaneously breaks the 
U(l) symmetry, which is accompanied by the appear- 
ance of an NG boson. This symmetry breaking pattern 
is summarized as follows: 

C/(1)b X 5[/(2)v (7(1)0 X (7(1) 
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For Nf = 3, we consider the 2+1 flavor case where the 
quark masses m„ = m^i ^ rus and the chemical potentials 
/^« 7^ 7^ Ms ill general. In this case, all flavor symme- 
try is broken explicitly leaving only U(l)xU(l)xU(l) in- 
variance which represents the conservation of each flavor. 
Once (e.g.,) d'j^s meson condensate is formed in addition 
to the nonzero chiral condensates, the symmetry is spon- 
taneously broken down to U(l) xU(l)ds, where U(l)ds 
leaves the condensate (d'j^s) invariant. We also find one 
NG mode corresponding to the spontaneous breaking of 
U(l) 0. The symmetry breaking pattern is 



C/(1)b X SU{3)y 



> U{1) X U{1) X Uil) 



K CO lid 



> C/(l) X (7(l)d.. 



(3) 



III. RANDOM MATRIX MODEL 



The chiral random matrix model[2^|23| is constructed 
based on the idea that the spontaneous breaking of the 
chiral symmetry is dominated by the low-lying Dirac 
eigenmodes, as is manifest in the Banks-Casher relation 
between the chiral condensate and the spectral density 
near zeroflJ]. The Dirac operator D is then truncated to 
be a matrix within the restricted space spanned by the 
(quasi-) zero modes, and the matrix elements are treated 
as random variables. Explicitly, in the representation 
where 75 = diag(l, —1), we have 



D 



iR + C 
ii?t + 



(4) 



where i? is a rectangular complex random matrix and C is 
a non-random matrix responsible for the effects of tem- 
perature and chemical potentials. Note that the block 
structure follows from the chiral symmetry {75, D} = 0. 

When the matrix R is not square but rectangular, the 
Dirac operator D has exact zero eigenvalues with \v\ 
being the difference between the numbers of the rows 
and the columns of R. This is the realization of the in- 
dex theorem for the topological charge and the number 
of exact zero modes in ChRM model. Hence, in order 
to include the anomaly effects, we need to deal with the 
non-square matrix properly. To this end, we categorize 
the zero modes into two species, the near-zero modes and 
the topological zero modes [l^. The former are assumed 
to be N left- and N right-handed low-lying modes gener- 
ated in gluon dynamics, while the latter are interpreted 
as the modes each of which is localized near one of 
instantons or iV_ anti-instantons in a gauge field con- 
figuration. Then the size of the matrix R is taken as 
(N + N+) X (N + N^) and the topological charge of the 



^ It is known that the numbers of the broken generators and the 
NG bosons can differ in some cases, e.g ., feromagnetism, kaon 
condensed phase with /ij = 0, etc. [28l. |29| . 



configuration is = — 7V_. We take the number of 
the near-zero modes 2N and the mean number of the 
topological zero modes {N±) to be proportional to the 
four volume V of the system. 

Concerning the matrix C which represents the medium 
effects, we introduce the effective temperature T and 
chemical potential fi f for a flavor / in the near-zero mode 
sector, while they are set to zero in the topological zero 
mode sector II9II23I: 



Cf = 



'(^/ +iT)ljv/2 









0^ 



ifif ~ ir)i^/2 I . (5) 

0; 



In the near-zero mode sector, fj,f±iT may be interpreted 
as the contributions from the two lowest Matsubara fre- 
quencies. The Dirac operator D with nonzero chemical 
potential 7^ is no longer anti-hcrmitian, though the 
partition function is invariant under /ij ^ ~M/- We 
stress here the fact that the absence of the medium effects 
in the topological zero mode sector resolves the unphys- 
ical suppression of the topological susce ptib ility of the 
original ChRM model (see discussion in |25l. [26|). This 
form of C/ may be understood physically as the fact that 
a topological zero mode localized near an (anti-) instan- 
ton is rather insensitive to the medium effects. 

Using the Dirac operator we deflne the partition 
function of the ChRM model with fixed N+ and 7V_ as 

r 

Zn^,n_ = dR c-^^''--^^' Y[ dctiDi^lf) + ruf) (6) 

/=i 

with a Gaussian weight for the matrix R. The typical 
scale of the chiral symmetry breaking is fixed by the pa- 
rameter S. The complete partition function is obtained 
by summing over the numbers of topological zero modes 
and iV_ with a distribution function P{N±) as 



Z 



RM 



E 



P(iV+)P(7V_)Zyv+,jv_ 



(7) 



The Poisson distribution would be appropriate for P{N±) 
in a dilute instanton system [13] ■ Adopting the Poisson 
distribution, however, one finds that the effective po- 
tential becomes unbounded from bclow|27|. Instead, we 
choose P{N±) to be a binomial distribution [Toj. 



P{N± 



N± 



p^±(l-p) 



'yN-N± 



(8) 



with parameters 7 and p. The binomial distribution mod- 
els the situation that, if we divide the four- volume V into 
7A^ cells, an instanton appears in one of the cells with the 
probability p, barring double occupancy. In other words, 
we introduced a repulsive interaction among instantons. 
This modification results in a bounded effective potential 
with a stable ground state[l9|. 
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After the standard bosonization manipulation, we can 
express the partition function at T = as 



Z^^\M, fi)= J dSe 



with the effective potential 

n{S; M, fi)^ _tr5^5_ -logdet 
7 



(9) 



^log|adet(S' + X) + l| 



S + M pL 

2 



(10) 



where the bosonic field S E C^i^^s is the order pa- 
rameter matrix corresponding to the bi-fermion field 



^fg 



and = diag(mu, md, . . . , mjvj) and 



fi = diag(/Ltu, (Id, ■ ■ ■ , fJ-Nf) are the mass and the chemical 
potential matrices, respectively. The term involving the 
parameters 7 and a = p/{l — p) breaks the U(1)a sym- 
metry S Se'-^ even in the A4 ^ limit [lil. [l7j . In the 
thermodynamic limit, N 00, the ground state is found 
as a solution of the saddle point equation, dfl/dSfg = 0. 
Using the solution S, the chiral condensate can be com- 
puted for a flavor / as 



d 



dm 



■n{S;M,fi) 



(11) 



By generalizing the mass matrix to the source matrix 
A4 = {s°- + ip'^)\°' / \/2 with the generators A" normalized 
tr(A°A'') ~ 2(5"'', other scalar (chiral) and pseudoscalar 
(meson) condensates are evaluated by differentiating in 
s° and p", respectively. 

In the preceding works [T9I [20j . we have studied the 
ground state of this model for Nf = 2 and 3 at finite 
temperature T and equal chemical potential /i. Regard- 
ing the anomaly parameters a and 7, wc have found that, 
for the large anomaly parameters, 0^7 ^ Y? , this model 
does not show chiral restoration, and therefore we should 
use the anomaly parameters in the region which allows 
the chiral phase transition. In the chiral limit the chiral 
phase transition at finite T and zero /i is found to be 
second (first) order for Nf =2 (3). This flavor- number 
dependence comes from the U(1)a anomaly term. Ex- 
tending to the fj, case, we find that the transition 
on the T-fi plane changes from the second order to the 
first order at a tri-critical point as // is increased in the 
Nf = 2 case. For Nf = 3 the first-order phase boundary 
separates the T-fj, plane into two regions. With increas- 
ing the quark mass, the thermal transition gets weak- 
ened and eventually turns to a smooth crossover, while 
the transition at larger fi than a critical point remains of 



first order. Then the T-fi phase diagram becomes simi- 
lar to the Nf=2 case with sma fl q uark masses. For more 
detailed discussions, see Refs. nj,M. 

IV. PHASE DIAGRAM: Nf = 2 CASE 



A. Effective potential 

We study the situation where /i„ 7^ fi^ with the de- 
generated quark masses m„ = = m. It is convenient 
to define the (averaged) quark chemical potential fiq and 
the isospin chemical potential /i/ as 



(12) 
(13) 



The order parameter matrix can be parametrized as 5* = 
Aa('/'a + ipa) with thc U(2) generators Aa (a = 0, 1, 2, 3). 
The usual ground state breaks the chiral symmetry spon- 
taneously, having nonzero chiral condensate ^q. At fi- 
nite ^7 the pion condensed phase where pi^2 are nonzero 
may be favored. Therefore, we adopt here the following 
Ansatz: 



P2 (pd 



(14) 



with real order parameters </>„, (j)d and pi,2- Notice that 
nonzero /i/ explicitly breaks the SU(2) isospin invariance 
down to U(l) invariance generated by A3, and then 0„ 7^ 
(j)d in general. Substituting this form into the effective 
potential ([T0|) . we obtain 

f^=^(0n + 0^ + 2|pp) 

- ^ log [(CTu + fJ'u){(Td - Pd) + 



log [{Ou - Aiu)('^d + Md) + \P? 



Uog[a{auCTd + \p\^) + l\ 



(15) 



where cr/ = <f>f -\- m. This model reduces to the one 
studied in Ref. 0] when the anomaly term is neglected, 
7a = 0, as it should. Note that the potential depends on 
pi and p2 only through \p\^ = p\-\- p2 due to the residual 
U(l) symmetry. Hereafter, we shall arbitrarily choose 
the meson condensate such that pi = and p2 ~ p, 
which breaks the U(l) invariance. Nonzero p ~ {{uj^d) — 
{d'y^u))/2 signals the pion condensation. 

Thc saddle point equations with respect to (pu ■, i>d and 
p respectively yield 



5 
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c^d + Pd 




7a(Td 
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f - ^J■d) - 




{(^u - Pu){<yd + Pd) - 


hp\ 
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. (o-„ - 


\- /Uu)(0-d - Md) - 
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1 
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_ (o-M H 
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{<Ju - Pu){o-d + Md) H 






- 1 



= 0, 
= 0, 
= 0. 



(16) 
(17) 
(18) 



Let us first classify tlie solutions with p = 0. Remember 
that when 07 = 0, the chiral condensates 0„ and (j)d 
decouple from each other. Thus, in the ideal case of to = 
0, we have four solutions by making combinations of 0y = 
and (f>f ^ 0, depending on the values of /iy. If aj ^ 0, 
however, the condensate 4'd cannot vanish exactly once 
4>u ^ because of the flavor mixing, while the trivial 
solution (j)u ~ (t>d ~ is still intact. Thus, in the chiral 
limit with 07 ^ 0, we expect four types of solutions, (i) 
ordinary phase with broken chiral symmetry, (ii) phase 
with restored chiral symmetry, (iii) & (iv) phases with 
nearly vanishing chiral condensate for only one of the 
flavors, u and d. We numerically find that the transitions 
between these four phases are of first order, and hence 
these four phases remain at small m 7^ 0. 

We next consider in what conditions the p 7^ solution 
appears. To this end, we expand the thermodynamic 
potential with respect to p at the ordinary ground state, 
dn/d(f>u,d ~ and p ~ 0: 

n{p;m,pf) = no{m,pf) + n2{m,pf)p'^-l . (19) 

Then one finds the coefficient up to 0{m,pj) as 

where (j>o is the scalar condensate at m = /z/ = 0. We 
see that the finite pion mass = toE^/^q disfavors the 
pion condensation, while the finite chemical potential pi 
favors it. From the condition D,2 = 0, the critical isospin 
chemical potential is found to be pj^ oc within this 
leading approximatior0. 

When m = i.e. = 0, infinitesimal pi selects out 
the ground state solution of ^„ = = and p ^ 
among the degenerate vacua. We find that from Eq. (fT8|) 
the non-trivial solution satisfies the relation 

= ^^uPd + f - ^^"j . (21) 



^ Note that wc use a different notation _A/,r for ttic pion mass in 
this model from the physical one m„ because the physical unit 
in the ChRM model is not fixed. 



For a'j = this recovers the solution of [8| 

pI = PuPd + ^ = - A^/ + ^ • (22) 

In the large p^ region, > solution disappears and the 
U(l) symmetry is restored. The second-order transition 
line is given by a hyperbola on the Pq-pi plane, obtained 
by setting p = in Eq. (|2T|) . This pion condensed region 
is enlarged by a factor (1 — a7/E^)~^ as compared to the 
case without anomaly. 

Before analyzing the ground state solution numerically, 
we remark here the symmetry of the effective poten- 
tial ([T0|) . First the effective potential is invariant under 
the respective charge conjugations, pu — > —pu and/or 
Pd — >■ —pd, when p — Q. In addition, the effective po- 
tential ()10p is invariant under u d. Reflecting these 
invariant operations, the phase diagram has an eightfold 
structure on the Pu-pd or Pq-pi plane. Nonzero con- 
densation p of the charged pion, however, breaks the re- 
spective charge conjugation symmetry leaving only the 
simultaneous one (pu, Pd) — >■ i—Pd,—Pu)- Together with 
the flavor symmetry, u -^r^ d, el fourfold structure remains 
in the phase diagram. On the Pq-pi plane, therefore, it 
is sufficient to investigate the phase diagram in the first 
quadrant, pq > and pi > 0. 



B. Numerical result 

In the top panel of Fig. [U wc show the phase diagram 
on the Pq-pi plane with parameters, a = 0.5, 7=1,1] = 
1 and TO = 0. The dotted lines denotes the first-order 
phase transitions with respect to the scalar condensates 
(j)u and 4>d when the possibility of the pion condensation 
is neglected. In this restricted case we find around the 
origin the ordinary phase where 4)u and 7^ 0. Allowing 
the pion condensation, we find that it completely covers 
the region of the ordinary chirally broken phase and it 
extends to the larger \pi \ region, as shown with the solid 
line in Fig.[T] Let us look at the small chemical potential 
region first. The ordinary phase with nonzero 0„ = (pd 
and /9 = and the pion condensed phases with p 7^ 
are coexisting along the line pi = \{pu — Pd) = in 
the chiral limit. With an infinitesimal pi, however, the 
ordinary ground state is totally rotated away to the pion 
condensed phase. This behavior was already found in the 
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chiral sigma model @. This is to be expected because the 
ChRM model is similar to the potential term of the chiral 
sigma model at small chemical potential /iy and small 
symmetry breaking mQ. Next, in the large-^y region 
in Fig. [Jl we see that the chiral symmetry restoration 
occurs. Increasing /i/ with fj.q kept small, we find a second 
order phase transition from the pion condensed phase 
to the chiral restored phase. On the other hand, if we 
increase /Zg along with the fixed ni line, we find a first- 
order phase transition from the pion condensed phase to 
the phase where one or both of the chiral condensates 
(pu and (pd melt away, depending on the size of fij. It is 
remarkable that the phase diagram with p = refiects 
the symmetry of {fiq,fii) {i^ii,fiq) with {4)m4>d) ^ 
{4>d,4'u), while in the real ground state with p 7^ 0, this 
symmetry no longer exists. As we mentioned, the phase 
diagram has the symmetry of (^q, /i/) — (±/Xg, i/i/). 

In the bottom panel of Fig. [T] we show the result with 
m = 0.1 7^ 0. Generally, finite m acts as an external 
alignment field for </>„, 4>d 0. So there is a competition 
between the two alignment fields m and iij. In a small 
fii region, the finite m wins and the ordinary phase with 
nonzero 1/)^ and (f>d appears pushing the pion condensed 
phase aside, as shown in Fig. [TJ If fij exceeds a critical 
value proportional to m^r, the pion condensation phase 
appears. It is estimated from Eq. (|20p with our model 
parameters as /x/c = \/(j)QlvQ/2 = 0.276 at = . We 
here remark the number of the phase transitions in small 
/i/ region where p = 0. When /ig is increased, the system 
experiences just one phase transition from the ordinary 
phase with 0„ > and > to the symmetric phase 
with (pu Q and ^ if |/z/| is sufficiently small. If 
l/x/l is large but still below the critical value, we have 
two transitions between the ordinary phase and the sym- 
metric phase; e.g., for /i/ > 0, the first one is from the 
phase of > 0, (pa > ^ to 0„ 0, <pd > 0, and the 
second is from 0„ ~ 0, (/>d > to 4>u 0, (pd ^ 0. In 
the study of the ChRM model without the anomaly ef- 
fect in @, a two-step phase transition is found. We find 
that, even with p = 0, the anomaly term mixes the <pu 
and <pd yielding a one-step transition. This mechanism 
is general and is found in an NJL model study fl3|. We 
expect that this is also the case in QCD. As a result, we 
have one triple point and two critical end points in the 
phase diagram shown in the bottom of Fig. [TJ (In the 
Nf = 3 case we will find more triple points and critical 
end points.) Varying the anomaly parameters a and 7, 
we have confirmed that the position of the triple point 
moves toward the larger /ig and p,i region as the strength 
of the flavor mixing is increased. 

Finally we comment on the large /i/ region regarding 
the possibility of BCS-like condensate. At large /i/ in 
QCD, large Fermi seas of (e.g.) u-quark and d-quark 
are formed, and the attraction between these quarks 



* A critical end point is the point where a critical line is truncated 
by meeting a first-order phase boundary. 




FIG. 1. Phase diagram for Nj = 2 on pq-fii plane in the 
chiral limit (top) and at m = 0.1 (bottom). Parameters are 
a = 0.5, 7=1, and E = 1. Nonzero condensates in the 
respective regions are indicated with the letters, p and pu.d- 
The first (second) order phase boundary is denoted in solid 
(dashed) lines. For comparison, the phase boundary in the 
case where the meson condensation is ignored, is denoted in 
dotted line showing the eightfold symmetry. 



may form a pion-like Cooper pairQ. However, since 
the ChRM model neglects the space-time dimensions, the 
physics of the Fermi surface does not exist and the su- 
perconducting phase is tricky. For possible extensions of 
the ChRM models to deal with the diquark condensates, 
see Ref. (soj . In our model without such an extension, we 
find a simple termination of the pion condensed phase at 
large \pi\. 
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V. PHASE DIAGRAM: iV/ = 3 CASE 

A. Effective potential 

We explore here the phase structure in the 2+1 flavor 
case with two mass parameters, m„ = = m and rrig, 
varying three chemical potentials fid and fis indepen- 
dently. They are recast to the quark chemical potential 
/ig, the isospin chemical potential pi and the hypercharge 
chemical potential fiy defined respectively as 



Pq= 2^tJ'U+ Pd), (23) 

= ^(Mu - A'd), (24) 
Mr = ^{Pu + Pd-^Hs). (25) 
I 



Note that fiy = Ps in our convention. When fii = 
fly = 0, three quark chemical potentials are equal to 
fiq. At nonzero fii and fiy, we allow the possibilities for 
the pion and kaon condensates in addition to the chiral 
condensates. We then apply the following Ansatz for the 
order parameter matrix, 

(0u Pud -Psu\ 
-Pud (pd Pds , (26) 
Psu -Pds 0s / 

where (pd, and 0s are the chiral condensates, p^d is 
the pion condensate and psu and pds arc the kaon conden- 
sates. With this Ansatz, we obtain the effective potential 
as an function of six order parameters, 

r! = fio + r!a, (27) 

where Q,q is the potential of the conventional ChRM 
modelQ, 



^^0 =\ {4>l + <fd + <fs+ 'ipls + 2pL + 2p,L) - \ ^og{{cjl - pl){al - fil){al ~ ij^) 

+ [Pdsi'^l - pi) + '^PdsPlui'^uCrd - PuPd) + 2pL('^« - Pl){'^dcrs - PdPs) + (cyclic perm, of u,d,s)] }^ , (28) 



and fia is the anomaly par10, 

^ 2 
^^a ^--\0g [a {(Tuddas + auPds + <^dPlu + <^sPld) + 1] 

(29) 

Because is a function of p^^, p^^ and pl^, we always 
have a trivial solution Pfg = for the saddle point equa- 
tions. On the other hand, once chiral and/or meson con- 
densates become nonzero, they act as source terms for 
the other chiral condensates owing to the anomaly term 
ria, and therefore cry = is no longer a solution. 

We confirmed numerically in the ChRM model that 
two or more meson condensates do not appear in the 
ground state at the same time. If only one type of 
the meson condensate pfg is nonzero, the ilo part of 
the potential becomes a sum of two contributions, = 
^ll^{pfg,(j)f,4>g) + rip (0;i), where fl^^ is nothing but the 
potential of the two-flavor (fg) ChRM model without 
anomaly, and fig corresponds to the single-flavor {h) 
ChRM model. In spite of this flavor separation, the phase 
diagram of Nf = 3 is different from that of Nf = 2 be- 
cause the anomaly term introduces the coupling among 
three flavors and we have a competition between the pion 
{fg =ud) and the kaon {fg =ds or su) condensed phases 



for being the ground state. Below, we just assume that 
only one component of the meson condensates becomes 
nonzero in the meson condensed phase and leave its proof 
as an open issu^. 

Preceding the numerical results, let us summarize 
three key points for qualitative understanding of the 
phase structure. 

(I) analogy with the chiral sigma model — As is shown 
in Ref. Q , the conventional ChRM model for small chem- 
ical potentials and quark masses is equivalent to the zero 
momentum part of the chiral Lagrangian. Hence the 
phase structure of the ChRM model in the small chemi- 
cal potential region must have a similar structure as the 
chiral sigma model. It is known in the sigma model with 
Nf = 3 that a second-order phase transition occurs from 
the ordinary chirally broken phase to the pion condensed 
phase or the kaon condensed phase at certain finite pi 
and py. The critical chemical potential is roughly esti- 
mated as Pi ~ and py ^ niK for the pion and 
kaon condensations, respectively. There is a competition 
between the pion and kaon condensed phases for finite 
Pj and py, and the phase transition between the two is 
found to be first-order. We also remark that the state 
with two meson condensates having nonzero values can 
not be even a meta-stable state in this analysis. 



Consequences of the anomaly mixing between the chiral and di- 
quark condensates are studied in Ref. [3ll . 



We note that coexistence of the p-wave pion and kaon conden- 
sates in nuclear matter was studied previouslv l32f . 
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(II) chiral restoration — In contrast to the chiral sigma 
model, the ChRM model includes the chiral restoration 
dynamics, which will result in a new kind of competition 
of two meson condensed phases in the region at large 
chemical potential. To illustrate the situation, let us con- 
sider the case in the chiral limit, where fiy > t'-i > t'-q ~ ^ 
and the i^°-condcnsed phase pus 7^ with 7^ is cho- 
sen as the ground state. We now further increase /zy with 

and /ij- fixed. This kaon condensed phase will simply 
remain as the ground state in the chiral sigma model. At 
sufficiently large py, however, chiral or meson conden- 
sates involving the s-quark become disfavored and the 
chiral symmetry in the s-quark sector will presumably 
be restored. Then there is a possibility for the remaining 
u- and d-quarks to form the pion condensate pud 7^ 0. 
Wc thus find a competition between the kaon condensed 
phase {pus 7^ O7 0d 7^ 0) a-nd the pion condensed phase 
[Pud 7^ 0, 0s = 0) at large p^y ■ We will see shortly in this 
section that the pion-condensed phase is indeed favored 
at lar ge n y. This competition exists also in the NJL 
model[l2|, and we expect that this mechanism is common 
for the models exhibiting chiral symmetry restoration at 
large chemical potentials. 

(III) anomaly effect — Mixing of the order parameters 
due to the anomaly tends to unify chiral phase transitions 
into a single one as is seen in the case of Nf = 2. In the 
absence of the meson condensates, the anomaly term Jla 
for Nf — 3 reduces to \og{aauCTd<7s + 1), which couples 
three chiral condensates. When one of the chiral conden- 
sates melts away, the flavor mixing among the chiral con- 
densates becomes small. With finite meson condensates, 
however, the anomaly term induces mixing among the 
chiral and the meson condensates. For example, once the 
pion condensate is formed p.^d ^ with rotating (/)„ and 
(\)d away, the mixing term a^ddOs becomes small but the 
term <i>sp\d ™ -'^'-l- generates new flavor mixing. Be- 
cause of this mixing, the chiral and pion condensates tend 
to vary cooperatively as the chemical potentials change. 
Furthermore, the anomaly term makes the ground state 
with nonzero 0s and pud condensates more stable as the 
term appears with a minus sign in the effective potential, 
which results in the extension of the meson condensed 
phase compared to the case without anomaly. 
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FIG. 2. Phase diagram for Nj = 3 on [ii-py plane with 
= in the chiral limit (top) and with nonzero quark masses 
— m,d — 0.02, nia =0.1 (bottom). In the bottom panel, 
the lower (upper) vertical narrow area represents the phase 
with nonzero cpu, (jjd and 03 (0^ and 0d). Other parameters 
and notations are the same as in Fig. [1] 



B. Numerical result 



We calculated the phase diagram with model param- 
eters S = 1, a = 0.5 and 7 = 1 in the chiral limit 
ruu = nid = TTis = as well as at finite quark masses, 
TO„ ~ nid — 0.02 and rus = 0.1. We show the phase dia- 
gram on the pj-py plane with pq = 0, the Pq-pi plane 
with py = 0, and the pq-py plane with pi = 0. 



1. pi~py plane 

We first present the phase diagram on the pi-py plane 
with zero quark chemical potential pq = in the chiral 
limit (top) and at finite quark masses (bottom) in Fig. [21 
Note that the effective potential has the symmetry of 
py —py and, moreover, pi o —pi if we change u 
and d fiavors simultaneously. We then present the result 
only in the first quadrant. 

Let us first focus on the case in the chiral limit. In 
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the small chemical potential region, we find the pion and 
kaon condensed phases. These two phases are separated 
by the first-order phase transition line /i/ = fiy the pion 
condensed phase appears when /i/ > fiy and otherwise 
the kaon condenses. This is also found in the chiral sigma 
modelQ. When chemical potentials are increased, we 
find that the regions of two mcsonic phases are exchanged 
on the diagram; the pion condensed phase appears when 
fij < fly and otherwise the kaon condenses. As already 
explained, this peculiar behavior is triggered by the chiral 
restoration; these meson condensed phases at the larger 
chemical potentials are accompanied by the melting of 
the chiral condensate of the remaining flavor. When ^y 
is further increased at small m, the pion condensed phase 
continues indefinitely because pud becomes insensitive to 
fiy once (ps disappears. On the other hand, pud continu- 
ously vanishes at some point as fij increases, because pud 
which has the isospin charge is affected by pi, whereas 
nonzero chiral condensate (j)s survives irrespective of pj 
for small py because it has no isospin. 

Without the anomaly term the vertical straight line of 
the second-order transition boundary seen in the large 
py region in Fig. [2] would extend down to the pi axis. 
This is because the effective potential without anomaly 
Hq becomes the sum of two contributions, the 0s part 
and the part involving (j)d and p^d, as mentioned be- 
fore, unless the kaon condensates have finite values. By 
comparison, we see that the anomaly coupling between 
(j)s and Pud makes the pion condensed phase more stable 
and extended to larger pi region. 

When quark masses are set to nonzero, the phase di- 
agram receives two qualitative modifications. One is the 
appearance of the chiral condensed phase without any 
meson condensate at small chemical potentials. This is 
the same result as in the chiral sigma model. The other 
is the configuration change of the four meson condensed 
phases. We find that the pion condensed phase extends 
from zero py to large py , cutting the two kaon condensed 
phases apart. This is because the chiral condensate (f)s is 
more favored for mg > ruu = md and the kaon condensed 
phases shrink. In the extended pion condensed phase, we 
find a first-order phase transition in 0^, which creates a 
small gap in p^d, too. 
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FIG. 3. Phase diagram for TV/ = 3 on pq-pi plane with py = 
in the chiral limit (top) and with nonzero quark masses 
(bottom). In the bottom panel, the horizontal narrow area 
is the phase with nonzero (pu, <i>d and Parameters and 
notations are the same as in Fig. [21 



2. Pq-pi plane 



We next show the phase diagram on the Pq-pi plane 
with py = in Fig. [31 In this case the effective poten- 
tial is symmetric in pq -H- —pq and in pj O —pi with 
exchange of the u and d quarks. We present the phase 
diagram again only in the first quadrant. 

Let us focus on the case in the chiral limit (top panel 
of Fig. 131). With infinitesimal pi, the ground state be- 
comes the pion condensed phase. Unexpectedly, how- 
ever, we find that a kaon-condensed phase appears in the 
region where pu becomes largest among three chemical 
potentials. This phase was not predicted in the chiral 



Lagrangian analysis. In this phase, the u-quark chiral 
condensate melts away 0„ ^ 0, and the remaining d- and 
s-quarks form the kaon condensate, pds 7^ 0. Note that, 
because rriK — 0, infinitesimal difference between pd and 
Ps makes the kaon condensed phase more stable than the 
chirally broken phase without a meson condensate. 

The pion and kaon condensed phases vanish contin- 
uously at large pj, respectively, to the phases with (ps 
and without any condensate. The kaon condensed phase 
extends to a larger pi region than the pion condensed 
phase. This may be understood from the fact that the 
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pion condensate with isospin charge 1 is twice as sensitive 
to /X/ as the kaon condensate with isospin 1/2. 

With finite quark masses we observe that two new 
phases are added, as is seen in the bottom panel of 
Fig. [3l The ordinary chirally broken phase appears in 
the small /i/ and region due to finite m,r, while finite 
rriK changes lower chemical potential part of the kaon 
condensed region into the chiral condensed phase with 
nonzero (fid and 0^. These two chiral condensed phases 
show the second-order phase transitions to the meson 
condensed phases at larger /i/, which is consistent with 
the chiral Lagrangian analysis [7|. At small /x/, we find a 
single transition from the broken to the symmetric phase 
along the /i^ axis, owing to the anomaly term. 

It would be instructive to compare this phase diagram 
to that of Nf = 2 on the fiq-^i plane. Although the 
phase diagrams are drawn in the same chemical poten- 
tial space, inclusion of the third quark fiavor changes the 
phase diagram drastically resulting in a new phase with 
the kaon condensation. If we take a limit of nis -> oo in 
the Nf = 3 phase diagram, the s-quark should decouple 
and the s-quark and kaon condensates disappear. This 
makes the phase diagram reduce to that of the Nf = 2 
case. 



3. l-iq-l-iY plane 

Finally, we address the phase diagram on the Hq~HY 
plane with jij = 0. In Fig. 2] we present the phase dia- 
gram in the chiral limit (top panel) and with finite quark 
masses (bottom panel). The effective potential is un- 
changed under the simultaneous exchanges of fj,q -H- —fiq 
and fiY ^ — /xy. Hence we restrict the phase diagrams 
in fly > region. The formation of the kaon condensate 
breaks the SU(2) isospin symmetry spontaneously. We 
choose arbitrarily the i^o-condensed phase {pds 7^ 0) as 
the meson condensed ground state. Note that no pion 
condensed phase appears in this diagram since we have 
set \fii\ = < m^r. 

In the chiral limit, the chiral and kaon (pion as well) 
condensed phases are degenerated at fiq = fxy = 0. 
With infinitesimal fiy, the ground state becomes the 
kaon condensed phase, and therefore the isospin sym- 
metry is broken spontaneously. Assuming the alignment 
by infinitesimal negative /i/, the chiral condensed state 
is completely rotated into the the K'^ condensed state 
where 0^ = 0s = and pds 7^ but the u-quark chiral 
condensate remains finite 0„ 7^ 0. When fiq is increased 
with fly = 0, two order parameters, pds and disap- 
pear at the same time at a critical value. Without the 
anomaly effect this phase transition would occur in two 
steps. 

We note that there is a first-order phase boundary in- 
side of the the kaon condensed phase pds 7^ 0. Across 
this line, (j)u jumps from a nonzero value to zero and pds 
changes also its value discontinuously. At one end of this 
line, we find a triple point where the two kaon condensed 
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FIG. 4. Phase diagram for N j = 3 on Hq-pv plane with 
pi — Oin the chiral limit (top) and with nonzero quark masses 
(bottom). Parameters and notations are the same as in Fig. [21 



phases and the <f>s ^ phase coexist. At the other end, 
we find a critical end point where two kaon condensed 
phases and the symmetric phase meet. 

Let us consider the finite mass effect. Because of the 
finite rriK, there appears the ordinary chirally broken 
phase in the small chemical potential region. The phase 
transition to the kaon condensed phase is of second order. 
Along the pq axis, we find the phase transition from the 
chiral condensed phase to the symmetric phase occurs in 
one step. The threshold py for the two-step phase tran- 
sition is relatively small compared to the corresponding 
value of Pi found in the Pq~pi plane (see Fig. [3]). 

This can be understood as follows: note that without 
the flavor mixing, the first order transitions along the 
Pq axis would occur at different values of pq when the 
quark masses are different. The U(1)a anomaly mixing 
acts to bind them together, and there is a critical value 
for the anomaly strength above which the phase tran- 
sitions occur in one step along the pq axis. The larger 
mass difference between u, d-quarks and s-quark makes 
this binding of the chiral transition more fragile against 
the external field py. This is the reason why we have a 
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smaller threshold value of fiy compared to the threshold 
value of fij. 

With the mass parameters rriu = md = 0.02 and rris = 
0.1 used in Fig.|31 the threshold value is found to be ac = 
0.279 . . . when 7 = 1 is fixed. In Fig. HI a = 0.5 > ac and 
we see the one-step phase transition along the fj,q axis. 

VI. SUMMARY 

We have investigated the phase structure of the ChRM 
model with 2 and 3 flavors on the plane of the quark 
chemical potentials, including the effects of anomaly. Dif- 
ferent chemical potentials for different flavors are real- 
ized in general situations because of the quark mass dif- 
ferences and electric neutrality, and the chiral conden- 
sates have also different values there. Moreover, when 
the chemical potential difference, which is flavor nonsin- 
glet, becomes large and comparable to the pseudo-scalar 
meson masses, we have a phase transition to a meson con- 
densed phase. We obtained a complex structure of the 
phase diagram with various orderings at first sight. Nev- 
ertheless, we can understand the diagram qualitatively 
based on the following three observations. 

The first point is the similarity to the chiral Lagrangian 
analysis, which explains the phase structure well below 
the chiral restoration. In the chiral limit for Nf = 2, the 
chiral condensed phase and the pion condensed phase are 
degenerated at fj,j = 0, and an infinitesimal ^/ rotates the 
chiral condensates into the pion condensate. With finite 
quark masses the chiral condensed phase survives up to 
yUj ~ m7r/2, and we find a second-order phase transition 
to the pion condensed phase. In the Nf = 3 case, a 
similar behavior is also found for meson condensations, 
but in addition we found a new competition between the 
pion and kaon condensed phases, bounded by a first-order 
phase transition. 

The second point concerns the chiral restoration, 
whose effect is not included in the chiral Lagrangian anal- 
ysis. For Nf = 2 and 3, we found that the pion and kaon 
condensates disappear continuously at large ^/ and fiy, 
respectively, while the phase transitions along the axis 
are first-order. The chiral restoration effect explains the 
peculiar appearance of the pion and the kaon condensed 
phases in the phase diagram for Nf = 3. In the chiral 
Lagrangian analysis, the pion condensation is always fa- 
vored if /i/ > fiY- However when /z/ is further increased, 
the chiral restoration for u-quark should occur and the 
pion condensation melts away. Then the kaon condensa- 
tion will be formed if the difference between fid and fig is 
larger than a certain critical value characterized by rriK- 
This argument also holds in the situation where fij < 
with exchanging the pion and kaon condensates. This is 
a remarkable inversion of the pion and kaon condensa- 
tion at high hj or /xy. Note that the argument does not 



depend on the details of the ChRM model. The same in- 
version may be found in other models which incorporate 
the chiral restoration dynamics, and therefore one may 
expect this also in QCD. 

The third point is about the anomaly effect. The most 
important effect is the mixing of the condensates. For 
Nf = 2 with finite quark masses, we find a single phase 
transition from the chiral condensed phase to the sym- 
metric phase along the /x, axis. Without the flavor mix- 
ing, the u- and d-quarks would show the chiral transitions 
independently at finite But these two phase transi- 
tions coalesce into a single transition via the flavor mixing 
term due to the U(1)a anomaly. There is the threshold 
lij for two-step phase transition, which becomes large if 
the anomaly parameters are increased. 

In the case of Nf = 3, the anomaly effect induces the 
mixing among not only the chiral condensates, but also 
the meson condensates, which makes likely that the me- 
son condensate and the chiral condensate of the remain- 
ing flavor show the discontinuity simultaneously as the 
chemical potentials are varied. In the mean-field models 
without the anomaly effect, where the flavor mixing is 
missing, the chiral phase transitions can happen in multi 
steps for finite ^/ and/or /zy even when the quarks have 
an equal mass. We expect that the chiral condensates are 
strongly correlated in QCD because of the anomaly mix- 
ing as well as the dynamics beyond the mean field, and 
it is unlikely to have the multi-step chiral restoration in 
QCD when fiq is increased with fij and being not too 
large. For large ^/ and/or ^y, however, there remains a 
possibility for chiral restoration to occur in multi steps, 
as shown in our model study. 

The ChRM model, which discards the space-time dy- 
namics but keeps the chiral symmetry only, is certainly 
a simplified model for QCD and will provide at most 
qualitative features of the QCD phase diagram. We have 
studied the response of the condensates to the chemi- 
cal potentials using this ChRM model with the U(1)a 
anomaly, and found interesting interplay between the chi- 
ral and meson condensates as well as the importance of 
the U(1)a anomaly effects. It will be intriguing to con- 
firm and elaborate the findings of this study by employ- 
ing other dynamic and microscopic models for QCD and 
direct simulations. 
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